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We investigate the statistical mechanics of chirality and biaxiality in liquid crystals through a 
variety of theoretical approaches, including Monte Carlo simulations, lattice mean-field theory, and 
Landau theory. All of these calculations show that there is an important interaction between 
cholesteric twist and biaxial order: The twist acts as a field on the biaxial order, and conversely, the 
biaxial order increases the twist, i.e. reduces the pitch. We model the behavior of chiral biaxial liquid 
crystals as a function of temperature, and discuss how the predictions can be tested in experiments. 
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In liquid crystals, there is a close connection between 
chirality, asymmetry under reflection, and biaxiality, ori- 
cntational order in the plane perpendicular to the di- 
rector. In the 1970s, Priest and Lubensky [l[ recog- 
nized that a cholesteric liquid crystal must have some 
slight biaxial order because of the difference between the 
directions along and perpendicular to the helical axis. 
Brand and Pleiner Q showed theoretically that chirality 
can smear out the transition between uniaxial and biax- 
ial phases, and Kroin et al. 0] confirmed this smearing 
experimentally in lyotropic liquid crystals. Later, Har- 
ris, Kamien, and Lubensky [4| developed a microscopic 
model of molecules interacting through classical central- 
force interatomic potentials, and found that cholesteric 
twist can only form if there are at least short-range biax- 
ial correlations between molecules. By comparison, in a 
system with quantum dispersive interactions, cholesteric 
twist can form even without such correlations [f|. 

In recent years, there has been a resurgence of interest 
in biaxial liquid crystals — driven in part by experimental 
reports of the discovery of a biaxial nematic phase in ther- 
motropic liquid crystals [H, 0] , and in part by prospects 
for usingbiaxial liquid crystals for fast-switching display 
devices 8]. For that reason, it is now important to re- 
examine the interplay between chirality and biaxiality in 
liquid crystals. The key issue is: How is the cholesteric 
pitch affected by biaxial order — either by long-range bi- 
axial order or by short-range biaxial correlations? 

In this paper, we investigate this issue through three 
theoretical approaches: (1) Monte Carlo simulations 
of a lattice model for chiral molecules interacting via 
anisotropic van der Waals forces. (2) Mean-field the- 
ory for the same lattice model. (3) Landau theory based 
on symmetry- allowed couplings between twist and biax- 
ial order. Through all three approaches, we calculate 
the cholesteric twist as a function of molecular chirality, 
molecular biaxiality, and temperature. These calcula- 
tions show that chirality acts as an effective field on the 
biaxial order, which changes the second-order uniaxial- 
biaxial transition into a rapid but nonsingular evolution. 
Conversely, biaxial order enhances the cholesteric twist, 
i.e. reduces the pitch, so that the pitch greatly decreases 
in the low-temperature, highly biaxial state. The calcula- 
tions also allow us to reconsider the relationship between 



twist and short-range biaxial correlations. Based on this 
theoretical work, we discuss opportunities for experimen- 
tal studies of chiral biaxial liquid crystals. 

For the simulations, we need a model molecular struc- 
ture that can exhibit biaxial order with or without chi- 
rality. Inspired by van der Meer et al. [H , we consider a 
structure with two ellipsoids arranged rigidly in the shape 
of the letter H. Each ellipsoid represents an extended, 
anisotropic charge distribution within the molecule. If 
the two ellipsoids are parallel, as in Fig. 1(a), this is 
an achiral biaxial structure, with a biaxiality character- 
ized by the separation h. By contrast, if the ellipsoids 
are twisted about the central connector, as in Fig. 1(b), 
this is a chiral biaxial structure, with a chirality char- 
acterized by the twist angle x of each ellipsoid from the 
parallel configuration. The interaction between any two 
ellipsoids on neighboring molecules is the van der Waals 
dipole-induced-dipole interaction. Hence, the total inter- 
action between two molecules i and j is the sum of four 
pairwise interactions among the constituent ellipsoids, 
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where e.; Q is the orientation of ellipsoid a on molecule i, 
and ri a jp — |r J( g — Vi a \ is the center-to-center distance 
between two interacting ellipsoids ia and jf3. 

Initially, we perform Monte Carlo simulations of achi- 
ral biaxial molecules. We simulate a simple cubic lat- 
tice of size 16 x 16 x 16, with a molecule centered on 
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FIG. 1: (a) Achiral and (b) chiral biaxial molecular structures 
studied in this work, (c) Cholesteric phase of chiral molecules, 
showing the macroscopic twist. 
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FIG. 2: Monte Carlo simulation results: (a) Uniaxial and biaxial order parameters as functions of temperature T, for achiral 
biaxial molecules with ellipsoid separation h — 0.24. (b) Complete phase diagram for achiral biaxial molecules, in terms of h 
and T. (c) Uniaxial and biaxial order parameters as functions of T, for chiral biaxial molecules with h = 0.24 and molecular 
twist angle x — 0.08. (d) Boundary twist angle <!> as a function of T. For achiral (% = 0) molecules, $ is locked at it, indicating 
that the system is not twisted. For chiral (x = 0.08) molecules, $ is not a multiple of it, indicating that the system is twisted, 
and the cholesteric twist increases as T decreases. 



each lattice site. In each Monte Carlo step, a molecule 
is randomly selected and its orientation is changed, fol- 
lowing the standard Metropolis algorithm. The uniaxial 
order parameter S and biaxial order parameter V are 
calculated as described by Bates and Luckhurst • For 
small ellipsoid separation h, we find a first-order transi- 
tion from isotropic (/) to uniaxial nematic (N u ), followed 
by a second-order transition to biaxial nematic (AT;,) at 
lower temperature. The temperature range of the biaxial 
nematic phase increases with h, as expected for board- 
shaped molecules. Figure 2(a) shows a sample plot of the 
order parameters for h — 0.24, and Fig. 2(b) shows the 
full simulated phase diagram. 

We now use the same approach to simulate chiral biax- 
ial molecules. In this system, we expect molecular chiral- 
ity to induce a cholesteric twist. This twist is generally 
not consistent with periodic boundary conditions. Hence, 
we use self-adjusting twisted boundary conditions in the 
z-direction, following the method of Mcmmer ll|. In this 
method, the boundary twist angle $ from the top to bot- 
tom of the cell is a free simulation variable, determined 



by the Monte Carlo process. In an untwisted system, $ 
must be a multiple of tt. Hence, the deviation of $ from a 
multiple of ir is a measure of the twist across the system, 
i.e. the inverse pitch. With this method, the simulation 
forms a cholesteric phase over a wide temperature range. 
A sample configuration showing the molecular orienta- 
tions along the z-axis is shown in Fig. 1(c). 

Using these simulations, we determine the uniaxial and 
biaxial order parameters for systems of chiral molecules. 
Figure 2(c) shows S and V as functions of temperature 
T for ellipsoid separation h = 0.24 and molecular twist 
angle \ = 0.08. At T — 3.4 there is a first-order transi- 
tion from isotropic to cholesteric, as seen from the jump 
in S. In the cholesteric phase there is a slight nonzero 
value of V, as expected from Ref. [l|]. As T decreases fur- 
ther, V gradually increases toward its maximum value of 
1. There is no phase transition between uniaxial and bi- 
axial, but only a nonsingular increase in V. Apparently 
the chirality acts as an effective field on the biaxial order, 
which smears out the N u -Nb transition. 

We also determine the boundary twist angle <& as 
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a function of T, as shown in Fig. 2(d). For achiral 
molecules, the boundary twist angle is locked at $ = 
7r, indicating that the system is in a uniform nematic 
phase, either uniaxial or biaxial. By contrast, for chiral 
molecules with x = 0.08, the results for $ show a twisted 
cholesteric phase. The cholesteric twist is substantial just 
below the isotropic-cholesteric transition, although the 
biaxial order parameter is very small there. The twist 
increases further as T decreases, and becomes largest in 
the temperature range that would be the biaxial phase 
for an achiral system. Thus, we see that the cholesteric 
twist and the biaxial order increase together, reinforcing 
each other, as T decreases. 

In addition to the long-range biaxial order, we mea- 
sure the short-range biaxial correlations as a function of 
distance between nearby lattice sites. Through most of 
the cholesteric temperature range, these biaxial correla- 
tions are very small; they do not become noticeable until 
slightly above the achiral N u -Nj, transition temperature. 

To compare with the simulations, we construct a 
Maier-Saupe-type mean-field theory for the same lattice 
model, following a method similar to our calculation for 
the flexoelectric effect [12] . Here, we assume the sys- 
tem has perfect order of the long axes of the molecules, 
but variable biaxial order and variable cholesteric twist. 
Suppose that site i has its long axis along the x-direction, 
as do the four neighbors in the xy-plane, while the two 
neighbors in the z-direction have long axes twisted about 
the z-axis. The long axis at site i is hi = (1,0,0), while 
the long axes of the neighbors are h± x = (1, 0, 0), h± y = 
(1,0,0), and h± z = R Z (±A6)(1, 0, 0), where A9 is the 
cholesteric twist from one layer to the next and R Z (A8) 
is the rotation operator about the z-axis. At each site, 
the molecular short axis must be in the plane perpendicu- 
lar to the long axis. Hence, the short axis at site i is b\ = 
(0, sin <f>i , cos (j>i ) , while the short axes of the neighbors are 
b± x = (O,sin0 ±x ,cos0 ±a; ), b± y = (0, sin <p± y , cos <f>± y ), 
and b± z = R Z (±A9)(0, sin cj)± z , cos 4>±z) , where the local 
angle (f> represents the azimuthal angle of the short axis. 

We now construct a distribution function for the local 
azimuthal angle cj>, which can be written as 



p{<£) = 



exp(C cos 2<fi) 
d(j> exp(C cos 2(j>) 



(2) 



where C is a variational parameter representing the ef- 
fective biaxial potential. It is related to the biaxial order 
parameter by V = cos(20)p(0)c?0. With this distri- 
bution function, the mean-field free energy per site can 
be written as the sum of energetic and entropic terms, 



F = (H) +k B T (log p), 



(3) 



where (H) is the average interaction energy of Eq. (p} 
between site i and its six neighbors. This free energy 
depends on two variational parameters (effective biaxial 
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FIG. 3: Theoretical results for cholesteric twist as a function 
of temperature (in units of interaction strength A): Monte 
Carlo simulations, mean-field theory, and Landau theory. 



potential C and cholesteric twist AO), two molecular pa- 
rameters (ellipsoid separation h and molecular chirality 
%), and temperature T (in units of interaction strength 
A). We numerically minimize the free energy over the 
variational parameters for each set of molecular parame- 
ters and temperature, to find the biaxial order and twist. 

The numerical mean-field results are consistent with 
the Monte Carlo simulations. For achiral biaxial 
molecules, the system has a uniaxial phase with V = 
at high temperature. As a critical temperature, it un- 
dergoes a second-order transition to a biaxial phase with 
f / 0, and the biaxial order parameter increases as a 
power law as T decreases. This achiral system is un- 
twisted, with A9 = for all T. By contrast, for chi- 
ral biaxial molecules, the system has a high-temperature 
cholesteric phase with a small but nonzero value of V. 
As T decreases, V increases gradually, without any phase 
transition between uniaxial and biaxial. The cholesteric 
twist A9 is substantial at high temperature, even when V 
is small, and it increases further as T decreases, as shown 
by crosses in Fig. 3. Thus, as in the simulations, we see 
that the chirality acts as a field that induces biaxial order 
and smears out the N u -Nb transition, and conversely, the 
biaxial order increases the cholesteric twist. 

For further insight into the relationship between 
cholesteric twist and biaxial order, we construct a Lan- 
dau theory for a chiral biaxial liquid crystal. As in the 
mean-field theory above, we suppose the system has per- 
fect uniaxial order along the local axis h(r), but variable 
biaxial order. The biaxial order can be described by the 
tensor By = V(bibj — CjCj), where V is the magnitude 
of the order, and b(r) and c(r) are the two principal 
axes orthogonal to n(r). The free energy can then be 
expanded in Bij and in gradients of n(r), to obtain 



F = ~K(dinj)(diTij) - Kq €ijknidjnk 
+|r(T - T UB )Tr(5 2 ) + isTr(B 4 ) 
-ue i1k Bnnidin k - wti^Bn^dhni 



T^Tr(B 6 ) 
(4) 
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In this expression, the first line is the Frank free energy 
for director gradients in a chiral liquid crystal, the second 
line is a power series expansion in , and the third line 
is a pair of chiral couplings between B t j and director 
gradients. If we now assume a cholesteric modulation of 
the form n = (cos qz, s'mqz, 0), b = (0,0,1), and c = 
(— singz, cos qz, 0), with an arbitrary twist wave vector 
q, the free energy simplifies to 
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In the limit of high temperature, where biaxial order is 
small and the s and t terms are negligible, we minimize 
this free energy over V and q to obtain 
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Kr(T - T UB ) ' 
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Equations © demonstrate that the twist q acts as a field 
on the biaxial order V, and conversely, the biaxial order 
increases the twist, and hence reduces the pitch. 

Instead of treating the Landau coefficients as purely 
phenomenological parameters, we can derive them from 
the lattice model presented in this paper, by expanding 
the free energy of Eq. ([3]) in powers of biaxial order and 
twist. Results of this calculation are shown by the solid 
line in Fig. 3. The predictions of Landau theory are con- 
sistent with simulation and mean-field results, except at 
low temperature where biaxial order is large and series 
expansion is unreliable. 

It is interesting to compare our results with Ref. 
which argued that short-range biaxial correlations are 
a key factor in determining cholesteric twist in systems 
with classical central- force interactions. We also find an 
important connection between biaxiality and cholesteric 
twist, but it differs from their argument in two ways: 

(a) Our model shows some twist even in the limit of 
no biaxiality. This result does not contradict Ref. [J], 
because our system does not have central-force interac- 
tions; it is consistent with Ref. [f| for quantum disper- 
sive interactions. However, it draws attention to the fact 
that most liquid crystals have quantum dispersive inter- 
actions, while central-force interactions are unusual. (To 
be sure, quantum dispersive interactions are derived from 
fluctuating microscopic central-force interactions among 
electrons, and these electrons might have some biaxial 
correlations. However, such correlations would be diffi- 
cult to observe in either experiments or simulations; nor- 
mal observations average over the fluctuations.) 

(b) Our results show there is not a cause-and-effect 
relationship between biaxiality and cholesteric twist; 
rather, there is a mutually reinforcing interaction be- 
tween them. Twist acts as a field on biaxial order, and 



conversely, biaxial order helps to increase twist. As tem- 
perature decreases, biaxiality and twist increase together. 



The theory presented here has implications for exper- 
iments on biaxial liquid crystals, either thermotropic or 
lyotropic. It should be possible to choose an achiral host 
that has a uniaxial-biaxial transition, and add a chiral 
dopant. The theory predicts that the dopant will induce 
a small twist (large pitch) in the uniaxial phase, but the 
twist will increase (pitch will decrease) as the uniaxial- 
biaxial transition is approached. For low dopant concen- 
tration the twist will diverge as (T — Tub) 1 , while for 
larger concentration the divergence will be more rounded. 
At the same time, the chiral dopant will smear out the 
uniaxial-biaxial transition. The chirality-induced round- 
ing of the uniaxial-biaxial transition has been observed 
in lyotropics Q, but has not yet been investigated in 
thermotropics. Moreover, to our knowledge, no experi- 
ments have yet examined the cholesteric twist around the 
uniaxial-biaxial transition in either thermotropics or ly- 
otropics. This should be a promising area for experimen- 
tal research, to further characterize the close relationship 
between chirality and biaxiality. 
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